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We have gained enough statistical precision to distinguish signal from noise in matrix elements of 
all operators relevant for the A7 = 1/2 rule in kaon decays and for the direct CP violation parameter 
e . We confirm significant enhancement of AI — 1/2 transitions observed in experiments, although 
a few large systematic uncertainties remain in our predictions: higher-order chiral corrections and 
lattice spacing dependence. The estimate of e'/e is further complicated by the problem of matching 
lattice and continuum operators. 



I. INTRODUCTION 

One of the poorly understood effects in low-energy phenomenology is the so-called AI =1/2 rule. Namely, kaon 
decays proceed at much higher rate through the AI =1/2 channel than through the AI = 3/2 one. In particular, for 
decays to two pions the following relationship between the amplitudes is experimentally observed: 

A(K -> (7r7r)j= ) _ A Q _ 

A{K -» (7r7r) 7=2 ) ~ A 2 { ' 

In the Standard Model, the electro-weak interactions and the short-distance part of strong interactions arc not 
enough to explain this predominance of AI = 1/2 transitions. The bulk of the ratio in Eq. [j] is attributed to the 
long-distance part of strong interactions. The long-distance effects are contained in the matrix elements (MEs) of the 
basis operators of the weak effective theory between hadron states. These MEs have remained mostly unknown due to 
their non-perturbative character. We compute them on the lattice with statistics more than sufficient to distinguish 
signal from noise for the AI =1/2 amplitude, thus allowing to check the prediction of the Standard Model and QCD 
against the experiment. 

In a related effort, we address the direct CP-violating parameter e', defined in 

<^+7r-|^ 2 ) 

(7r+7r-|Jfi) ' [ 1 

where K\ (K2) are CP-even (odd) strong interaction eigenstates. Our goal is to predict the value of e' in the Standard 
Model in order to see if it is different from the Superweak theory and to compare it with experiment. As for the exper- 
imental situation, Fermilab's KTEV group has recently announced their latest result of Re(e'/e) = (28 ± 4.1) • 10 -4 , 
which is consistent with the old result from CERN's NA48 experiment ((23 ± 7) • 10~ 4 ), although hardly consistent 
with the Fermilab's previous result of (7.4 ± 5.9) • 10 -4 . The e'/e value will be even more refined soon by including 
more data from Fermilab and by complementing it with results from two other high-precision experiments at CERN 
and 4> factory at Frascati. The theory is lagging behind these experiments, mostly due to the lack of knowledge of 
MEs containing long-distance dynamics. 
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II. SIMULATION DETAILS 

We work within the effective theory obtained from the Standard Model by integrating out the W boson and t, b 
and c quarks p|. The effective Hamiltonian is written in terms of linear combination of basis four-fermion operators: 

H$ = -^V ud V* s ]T + rvifji)] OM , (3) 

where Zi and yi are Wilson coefficients (currently known at two-loop order) and r = ~VtdV t */V u dV* s . In this work 
we seek to compute MEs (Tnr\Oi\K), which are necessary for estimation of both A^ and A2 amplitudes and e'. 
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TABLE I. Lattice parameters 



Nt 
j 


0= \ 


Lattice size 


L, fm 


Number of 


Quark masses 










configurations 


used 





6.0 


16 s x (32 x 4) 


1.6 


216 


0.01 — 0.05 





6.0 


32 3 x (64 x 2) 


3.2 


26 


0.01 — 0.05 





6.2 


24 3 x (48 x 4) 


1.7 


26 


0.005 — 0.03 


2 


5.7 


16 3 x (32 x 4) 


1.6 


83 


0.01 — 0.05 



For technical reasons, two-hadron states are very difficult to put on the lattice ||. We compute {ir\Oi\K) and 
(0\Oi\K) and use chiral perturbation theory to recover (7nr\Oi\K) as follows 0: 

1 9 

(vr+vr-IO.^ ) = (tt+IO, - a i O aub \K+) ■ JL (4) 

(Pt, ■ PK)J 

where 

(0\Oi\K°) ... 
ai (0\O sub \K°) ■ (5) 

With staggered fermions, there is only one (lower-dimensional) operator that needs to be subtracted non-perturbatively 
in the manner shown above, namely 

Osub = {md + m s )sd + (rrid - m s )s^5d . (6) 

It should be kept in mind that this procedure is based on the lowest-order relationships in the chiral perturbation 
theory, and thus may be significantly biased. In particular, it ignores the final state interactions between pions among 
other higher-order corrections. 
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In Table [j] we show the simulation parameters. We are using a quenched j3 = 6.0 ensemble, and compare it with a 
quenched (3 — 6.2 ensemble to study cutoff dependence, with a (3 = 5.7 ensemble with two dynamical quark flavors 
(with comparable cutoff size) to check the magnitude of quenching effects, and with another quenched [3 — 6.0 
ensemble to check the finite volume dependence. 

Shown in Fig. [I] are the diagrams we are computing. The "eight" diagrams, along with the two-point function and 
"subtraction" diagrams, are relatively cheap to compute. In contrast, the "eye" and "annihilation" diagrams have 
been known in the past to be noisy and therefore difficult to analyze. In the present study we have gained enough 
statistics to compute them with reliable accuracy. 

We use staggered fcrmions and gauge-invariant, tadpole-improved operators throughout the simulation. The masses 
of all light quarks are equal. Other details and exact expressions for various quantities can be found in our upcoming 
paper @. 

III. RESULTS OF THE SIMULATION: AI = 1/2 RULE 

The main results of this work are shown in Fig. |[ The dependence of isospin amplitude ratio on the kaon mass 
is quite dramatic. This is due to the behaviour of Re^b amplitude, shown in Fig. || (for dynamical ensemble). The 
Re^4o amplitude is rather weakly dependent on the kaon mass (Fig. [|). Dependence of Re^lo amplitude on the cutoff 
size is shown in Fig. |^. The finite volume dependence and the effect of quenching are found small compared to noise. 

In order to compare the theoretical prediction for the amplitude ratio with experiment, it is necessary to specify 
a mass point for extrapolation in Fig. || If higher-order chiral relationships were used to obtain these amplitudes, 
they would provide a natural mass point. At this stage, all we can do is suggest that the mass point is somewhere 
between kaon and pion mass. This would give us a result approximately consistent with experiment, especially for 
the dynamical ensemble. Independent of the exact mass point, a significant enhancement of AI = 1/2 over AI = 3/2 
transition is evident. 




FIG. 2. The ratio Re^/Re^lo versus the meson mass squared for quenched and dynamical ensembles. The dynamical 
ensemble data were used for the fit. The horizontal line shows the experimental value of 1/22. The vertical line corresponds to 
the physical kaon mass. The error bars show only the statistical errors. 
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FIG. 3. Re^2 for the dynamical ensemble. The horizontal line is the experimental value of 1.23 GeV. 
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FIG. 4. Re^o for quenched ensembles plotted against the meson mass squared. The upper group of points is for ensembles 
>i and Q2, while the lower group is for Q3. Only statistical errors are shown. 
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FIG. 5. ReAo for quenched ensembles plotted against lattice spacing squared. The horizontal line shows the experimental 
result of 27.2 ■ 10 -8 GeV. Only statistical errors are shown. 



IV. DIFFICULTIES WITH e' 

In addition to the difficulties mentioned above, estimation of e' is further complicated by the issue of matching 
lattice and continuum operators. Usually such matching is done using lattice perturbation theory. For some oper- 
ators important for e' , the perturbation theory is unreliable. This is a serious obstacle for estimating e' , especially 
considering that it is an extremely fragile quantity, dependent on several subtle cancellations. What is needed is a 
non-perturbative matching procedure, similar to the one described in Ref. Q. 

As a temporary step, we have adopted a partially non-perturbative matching procedure, based on computing bilinear 
renormalization constants . We can obtain only an approximate guide to the size of four-fermion operator matching 
coefficients, since a few ad hoc assumptions are made. If this procedure is true, we get a surprising result: e'/e < 0. 
This result and the validity of our assumptions should be verified by full non-perturbative operator matching. 



V. CONCLUSION 

We have computed MEs of all basis operators appropriate for predicting the magnitude of AI = 1/2 amplitude 
enhancement and s' /e, on a number of gauge ensembles with reasonable statistical accuracy. A few systematic 
uncertainties still preclude a definite result for either of the two quantities. For AI = 1/2 rule, the main uncertainty 
is due to unknown higher-order chiral corrections, while for e'/e there is an additional uncertainty due to the difficulties 
with operator matching. We observe that within this sizeable uncertainty, the prediction for Re^lo/Re^ ratio is of 
the same order as the experiment. 

The computations were done at the Ohio Supercomputing Center and NERSC. We thank Columbia group for access 
to the dynamical configurations. 
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